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[57] ABSTRACT 

The present invention provides for sample timing adjust- 
ment using channel impulse response which is particularly 
weU adapted for use in a receiver for digital cellular mobile 
telephones. The sample timing adjustment includes storing 
samples in time slots, generating estimates of the transmis- 
sion channel impulse response, and using a subset of the 
samples received during the time slot, and generating error 
measurements of the degree to which the estimated channel 
impulse responses match the actual channel impulse 
responses. The different measures use simultaneously 
recorded samples having different time offsets. At least one 
sample used for each measure is advanced in time and one 
sample used for each measure is retarded in time relative to 
the samples used for generating the first measure. Then, a . 
sample time setting that minimizes the error measurement is 
determined and the sampling is adjusted to reflect the newly 
determined sample time setting. Time-reversed operation 
further enhances performance and permits implementation 
with reasonable complexity for a standards compliant 
mobile receiver. 

6 Claims, 6 Drawing Sheets 
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RECEIVER SAMPLE TIMING ADJUSTMENT 
USING CHANNEL IMPULSE RESPONSE 

This application is a division of U.S. patent application 
Sen No. 07/722,440, filed Jun. 27. 1991, U.S. Pat. No. 
5,263,026, assigned to the same assignee as the present 
invention. 
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BACKGROUND 

. The present invention relates' generally to digital cellular 
communications, and more particularly, to a maximum 

L likelihood sequence estimation based equalization method 
for use in mobile digital celluiar receiv ers." 

^Communication channels in the cellular chvifbnment 
|Cominonly impose a combination of distorting effects on 
(^transmitted signals. Rayleigh fading, where a signal's per- 
ceived power level rises and falls rapidly over a wide range, 
results from the combination of signals that have traversed 
paths differing in length by at least a significant fraction of 20 
a wavelength (i.e., about 30 cm. for cellular). Differences in 
path transmission times that approach the time taken to 
transmit a symbol result in a second problem called delay 
spread. 

Delay spread results in reception of multiple delayed 25 
replicas of a transmitted signal. Each Rayleigh faded replica 
has randomly distributed amplimde and phase, and the rate 
at which this complex quantity varies is constrained by the 
Doppler bandwidUi associated with a vehicle's speed. In a 
frequency nonselective environment, the sampled outputs of 30 
a receiver's matched filter^provides uncorrelaied estimates 
of the transmitted data. ^As~such^^in^tcrms^of. discfiete'time \ 
. samples, the channel has exhibited an impulse response 
proportional to a delta function . With delay spread, on th e 
other hand, the discrete time ch"ahne] impulse^esponse is 
extended to introduce energy at a number of symbol times . 
c The^ffect of the channcl-on the transnaitted signal, in turn , 
may be viewed as'the mnvoluiig|i of the transmitted infor m 
ination~with~thc cSaSS T^mpuI^'Icsponse ^^ The' channel;^ 



Ui^ efbre^^^^ s_^^axg^gU^^^^fi^g^^ '~ 40 

TmHeadls to the possibility of estimating the transmitte d 
information through the use of methods analogous to typical 
decoding of convolutional codes , lc., maximum luceiinood 
sequence estimation technique s. Unlike the more widely 
applied forward error correction decoding environment, the 45 
details of the encoding process in a leveree enor correction 
decoding environment, are not known a priori by th e 
receiver. Issues related to the need to estimate the form of "S e 
encoding process arc addressed by this inventi on. 

For the North American digital cellular system, a number 
of documents define the standards of implemented compo- 
nents; With respect to this invention, the following arc of 
interest: "Dual-Mode Mobile Station — ^Base Station Com- 
patibility Standard" denoted here as IS-54, [EIA/TIA Project 
Number 2398. Rev. A, January 1991] and "Recommended 
Minimum Performance Standards for 800 MHz Dual-Mode 
Mobile Stations", denoted here as IS-55, [EIAATA Project 
Number 2216, April 1991]. 

It would therefore be desirable to provide an enhancement ^ 
to the processing performed by equalizers for use in mobile 
telephones that provides for system complexity reduction 
and that provides for better performance in a fading channel. 



cellular mobile telephones. A number of distinct aspects of 
the design are considered novel, including: non-real time 
operation mode that permits and exploits **time-reversed" 
equalization, which significantly er^ances bit error rate 
performance ; the use of maxirnum likelihood sequence 
estimati on; use of variable coefficient least mean square 
Sacking dimng estimation of the transmission channel 's 
impulse respon se; and integrated symbol timing adjustment 
and carrier tracking algorithms. 

More particularly, die present invention comprises a 
method of processing received symbols including known 
preamble data and transmitted data, which method compen- 
sates for the effects of a power fade caused by a frequency 
selective fading channel. One specific aspect of the present 
invention comprises a method of processing samples 
received from a delay-spread, fading chaimel, which 
samples are associated with a block of data transmitted 
within a time slot, and wherein the method is adapted to 
make data decisions, and comprises the following steps: (1) 
Storing samples received during the time slot; (2) Estimating 
the location within the time slot at which decision errors are 
most probable, by determining the location in time of the 
maximum fade depth in the transmission charmel impulse 
response; (3) Processing the stored samples, starting with 
the first received sample and proceeding beyond the location 
of the maximum fade depth, using a predetermined maxi- 
mum likeUhood sequence estimation procedure to generate 
estimates of the transmitted data; (4) Processing the stored 
samples, starting with the final received sample and pro- 
ceeding in a reverse direction with respect to the time 
sequence in which the samples were stored, beyond the 
above-determined location, using the maximum likelihood 
sequence estimation process to generate cstunates of the 
transmitted d at a; (S) /Simultanenus with 'the" above" two? 
ynfocessing step^t ^enCTaUng Rstimatfts o tmecna ractensDc^ 
l ot the tr^ ^^^ ^ chaimel impulse res^mse^mcn are 
^ea m the maximum like'EEOTH^eciuence esumaaon nro ^ 
^S3: | f6) Processing the output of the preceding estiniate 
generating steps to generate data decisions. 
The step of generating the channel impulse response — 



SUMMARY OF THE INVENTION 

The present invention provides for an equalization 
method that forms the nucleus of a receiver for digital 
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estmiates typically comprises using variable tap coefficients 
that are determined by estimating tap settings within the 
estimated charmel impulse response by minimizinE the 
square of the difference between actual received samples 
and those synthesized by passing known transmitted sienaj s 
flirough the estimated channel. The processing is done in an ' 
ita:auve manner by combining previous estimates of chan - 
nel impulse response and new estimates thereof based on 
recent information, and by varyinfl the ratio of the contri - 
l>utiops from the previous and new estimates as a frinction of 
location within me time sioL " 

The method may further comprise a symbol timing adjust- 
ment procedure comprising the following steps: (1) Using a 
subset of the samples received during a time slot, generating 
an error measurement comprising a measure of the degree to 
which the estimated channel impulse response matches the 
acmal chaimel impulse response; (2) Generating a plurality 
of similar measures utilizing simultaneously recorded 
samples having different time offsets wherein at least one 
sample is advanced and one sample is retarded in time 
relative to the above measure; (3) Searching for a bit time 
setting that minimizes the above error measurement and 
adjusting the sampling to reflect the newly determined bit 
time setting. 

The method may also further comprise a carrier offset 
tracking method comprising the following steps: (1) Record- 
ing at least two samples of at least one tap within the 
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estimated channel impulse response at selected symbol 
locations within the time slot during the equalization pro- 
cess; (2) Generating frequency offset estimates at each of a 
plurality of time slots by observing the phase difference in 
each time slot between the at least two samples; (3) Com- s 
billing this plurality of frequency offset estimates by using a 
filtering process to generate a precise frequency offset 
estimate; (4) Adjusting a conurollable frequency source to 
compensate for the precise frequency offset estimate. 

Also, the present invention also uses a forward maximum 
likelihood estimation method of processing samples 
received from a delay-spread, fading channel, which is 
adapted to make data decisions. This method comprises the 
following steps: Storing samples received during the time 
slot; Processing the stored samples, starting with the first 13 
received sample and proceeding beyond the last received 
sample of the time slot using a predetermined maximum 
likelihood sequence estimation procedure to generate esti- 
mates of the transmitted data; Simultaneous with the above 
processing, generatinj^ estimates of the characteristics of th e ^ 
transmission channel impulse response which arc used in 
maximum likelihood sequence estimation process; Process- 
ing the outputs of the preceding steps to generate data 
decisions. 

The maximum likelihood sequence estimation method 
offers significant performance advantages when compared 
with altemative equalization options, such as decision feed- 
back equalization, for example. With respect to meeting 
industry defined standards for performance of digital cellular 
telephones used in fading envlrorunents, the performance of 
the equalizer incorporating the maximum likelihood 
sequence estimation process of the present invention stands 
out as superior. It is tiie only known process that is capable 
of meeting, or even approaching, the current digital cellular 
mobile telephone specifications. Time reversed operation 
further enhances p^onnance and permits implementation 
with reasonable complexity of a standards compliant mobile 
receiver. 

y 

40 

BREBj^^^RIPTION OF THE DRAWINGS 

The various*fcatures and advantages of the present inven- 
tion may be more readily understood with reference to the 
following detailed description taken in conjunction with tiie 45 
accompanying drawings, wherein like reference numerals . 
designate like structural elements, and in which: 

nG. 1 illustrates the problem associated with reception in 
a mobile environment with a fading channel and one aspect 
of the solution provided by the present invention; 

FIG. 2 is a block diagram of a digital cellular mobile 
telephone receiver incorporating a maximum likelihood 
sequence estimation based equalizer in accordance with the 
principles of the present invention; 

FIG. 3 shows the processing performed in die maximum 
likelihood sequence estimation based equalizer of FIG. 2; 

FIG. 4 illustrates a graph of relative error versus bit tinting 
offset that illustrates one problem overcome using the maxi- 
mum likelihood sequence estimation based equalizer of the 
present invention; 

FIG 5 is a flow diagram illustrating the processing 
performed by the equalizer of the present invention to 
implement carrier frequency offset compensation; and 

FIGS. 6A and 6B are flow diagrams illustrating the 
processing performed by the equalizer of the present inven- 
tion to implement bit timing. 
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DETAILED DESCRIPTION 

Referring to tiie drawing figures, FIG. 1 illustrates the 
problem associated with reception in a mobile environment 
having a fading channel and one aspect of the solution 
provided by the present invention. FIG. 1 shows a graph 
showing received power level from a typical fading channel 
versus time. The location of the power fade is shown relative 
to a typical time sloL The time slot is shown enlarged and 
includes a preamble and a coded digital verification color 
code field (CDVCC), which comprise known dam that is 
used to initialize a receiver system employing the equalizer 
of the present invention. At the lowest portion of FIG. 1, 
equalization processing in accordance with the present 
invention is illustrated with die arrows "A" and "B*\ in 
which tiie equalizer of the present invention processes 
forward and time reversed computations through the loca- 
tion of the power fade in order to accomplish the objectives 
of the present invention. This will be more fully described 
below with reference to FIGS. 2 and 3. 

FIG. 2 is a block diagram of a digital cellular mobile 
telephone receiver system 20 incorporating a maximum 
likelihood sequence estimation based equalizer 21 in accor- 
dance with the principles of the present invention. The 
system 20 comprises an amplifier 22 whose output is 
coupled by way of a downconverter, comprising a frequency 
source 23 and a mixer 24, to an analog filter 25. An analog 
to digital converter 26 is coupled to the analog filter 25 in 
order to digitize the downconveited data. A matched filter 27 
is coupled between the analog to digital converter 26 and the 
equalizer 21 of the present invention. The equalizer 21 
comprises a memory 30, a 4-state equalization trellis 31 that 
is adapted to calculate maximum likelihood sequence esti- 
mation metrics, a channel iimiulse response estimator 32, 
and an equalizer control circuit 15 ! 

A serially coupled AGC circuit 35 and gain control circuit 
38 are coupled to the amplifier 22. The equalizer control 
circuit 33 is coupled to an output of the matched filters 27 
and is coupled to an input to the frequency source 23. 
Symbol sampling (bit timing) time control circuitry 37 is 
coupled to the equalizer control circuit 33 and the acquisi- 
tion circuit 36 and provides control signals to die analog to 
digital converter 26. The output of the matched filters 27 is 
coupled to the AGC circuit 35 and the acquisition circuit 36 
and to the equalizer control circuit 33 that is employed to 
control the frequency source 23 and provide traiiung data for 
use in irutializing the equalizer 21. 

In operation, a partially filtered IF sip ial with a renter 
frequency of 85.05 MHz enters the gain controllable ampli - 
fier 22. The resulting signal is then downconverted using the 
frequency source 23 and the mixer 24 to 461.7 kHz, Hiis 
signal is then filtered using a nanow analog filter 25 to reject 
most of the received signals outside die 30 kHz band of 
interest The resulting signal is then sampled and converted 
to 8-bit digital samples using die analog to digital (A/D) 
converter 26. A 16 tap fractionally spaced digital FIR filter 
27 then performs matched filtering to produce symbol 
spaced samples which enter the equalizer 21. Temporally 
offset matched filters 34 that are substantially the same as the 
matched filters 27 are provided for use by the symbol timing 
control circuit 37, via the equalizer control circuit 33. 

The principles of maximum like^ il^ ^d sequence estima - 
tion employed in the equalizer 21 have been described in 
teclmical literature starting in the early I970's. A usef ul 
ouUine is found in "Adaptive Maximum-Likelihood 
Receiver for Carrier-Modulated Data Transmission Sys- 
tems", by G. Ungerboeck, IEEE Trans, on Communications, 
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Vol. COM-22, pp. 624-636, May 1974. Another description are useful for training. Assuming that a channel impulse 
of the maximum likelihood sequence estimation technique is response tracking breakdown occurs, this approach mini- 
provided in the reference "Digital Communications —2nd mizes the number of affected symbols by predicting the 
Edition.", by J. G. Proakis, 1989, pp. 610-642. failure point and avoiding equalization beyond that point 
The maximum likelihood sequence estimation process is 5 At 100 km/hr, which is the maximum speed specified in 
outlined as follows. The channel has an impulse response IS-55, which describes the mobile unit minimum perfor- 
containing significant energy in, say, N symbols. Assum e mance requirements, the average time between fades are on 
that the transmitter sends a sequence or symbols, much the order of 12 milliseconds. Given time slot durations of 
longer than N. The transmitted sequence may be described. about 6.7 milliseconds, there is only a small possibility of 
as the transitions between states, where each state coire^ iq two significant fades occurring within a time slot. However, 
sponds to a group of N-1 transmitted svmbols. The states! very close to the center of the slot is the coded digital 
therefore, correspond to overlapping groups of transmitted verification color code field. Even after a channel impulse 
symbols. In consecutive states, therefore, all but one con- response tracking breakdown, the channel impulse response 
stituent symbol are the same, and the possible transitions estimator 32 is very likely to recover during processing of 
between states are coirespondinfily constrained . As each „ the coded digital verification color codes due to the certainty 
sample is received, the equalization trellis 31 considers of the , transmitted data. Hence, the underlying period for 
every possible sequence of N svmbols that could hav e" which multiple fades are a concern is around 3.5 milUsec- 
contnbutcd to itTvalue. bv convol^ng that sequence witi i- on^s. The chance of more than one deep fade occumng 
. the estimated channel impulse respons e; For each hypoth- during this time is very low. ConsequenUy, tmie-reversed 
r esized sequence, the result of the convolution corresponds , equalization improves bit error rate performance m the 
, or fails to correspond, in some wav Cdefined bv a sUtistiC digital cellular envu*onment. 

called a metric) to the measured sample . On an individual The present equalizer 21 uses a 4-state architecture, 

basis, the hypodiesizcd sequence with the closes^ match to corresponding to N=2, where N is the length of the estimated 

Uie measured sample (the best metric) is the most likelv to channel impulse response. This choice assumes that the 

have been transmitted . However, over many samples and 25 energy in two (symbol-spaced) samples of the charmers 

under the c onstraint lhat only certain state transitions are impulse response dominates. To avoid channel impulse 

possible , t he path (sequence of states) with die minimum response tracking breakdown problems, reverse equalization 

cumulative metnc has maximum likelihood, and this is what is used for those symbols following the minimum power 

t he decoder selects^ point in a received time slot. 

The system 20 has no a priori knowledge of the form of 30 More specifically. FIG. 3 shows the processing performed 

the encode r employed in the transmitte r. Performance of the in the maximum likelihood sequence estimation based 

equalizer 21 therefore depends on the accuracy of the equalizer 21 of FIG. 2. The first step involves finding the 

estimate of the encoder's state , t he channel impulse response location of the power fade (box 51) in terms of symbol 

CCIR). FIG. 2 also shows the signals used in estimating the number. Processing starts in the forward direction toward 

cEannel'm^julse response.. The objective is to estimate the 35 the location of die power fade. The symbol number is set to 

form of the transversal finite impulse response filler that 0 (box 52), and then incremented (box S3). A decision is 

would take as its input the transmitted information symbols made whether the symbol then processed is a training 

{a(n)}, and produce at its output the samples taken from the symbol (box 54). If the symbol encountered is a training 

matched filter, {z(n)}. During the transmission of preambles symbol, then training data is inserted (box 57). If a training 

and coded digital verification color codes, the receiver 40 symbol is not processed, then the equalization trellis is 

knows the values of {a(n)}. However, at other times, only employed to generate metrics and, if possible, a decision 

the estimated values {aw(n)) are available for use in the (box 55). This is accomplished using equations outiined 

channel impulse response estimation process. T his depen- below. Then it is determined if a decision has been made 

dencc leads to a significant performance-degrading possi- (box 56). If a decision has been made, then an estimate of the 

bility. If decision errors emerge from the equalizer, and these 45 channel impulse . response is generated (box 58). If die 

are then used to update the estimate of the charmel impulse decision is not made, or once the channel impulse response 

response, tixen further decision errors become more probable estimate has been generated, then the symbol number is 

leading in a circular fashion to further decision errors and compared to the location of the power fade plus a predeter- 

breakdown of the equalization process. This phenomena is mined number of additional symbols Cox 59). Processing is 

referred to as a "chaimel impulse response tracking break- 50 then repeated by incrementing the symbol number (box 53) 

down'*. Such difiBculties are most likely to arise at die and repeating steps (boxes 54-59) until die fade location 

periods of minimum signal-to-noise ratio, or when the plus a predetermiried number of additional symbols has been 

received signal power is at its minimum during reception of reached 

a slot Once the desired symbol location is reached in (box 59), 

Widiin the IS -54 standard, which describes the interface 55 then processing is performed in the reverse direction starting 

between mobile and base eqidpment for North American widi the preamble of the next succeeding time slot, namely 

digital cellular systems, each information time slot is pre- symbol number 177, for example. The symbol number is set 

ceded by a known sequence, designated as the "preamble". to 178 (box 62), and then decremented (box 63). A decision 

As viewed by the receiver, therefore, information in the time is made whether the symbol then processed is a training 

slot is bounded on both sides by known sequences; the 60 symbol (box 64). If the symbol encountered is a training 

preamble for this slot and the preamble for the subsequent symbol, then training data is inserted (box 67). If a training 

slot. Consequently, this equalizer 21 is adapted to mitigate symbol is not processed, then the equ^ization trellis is 

the effects of a channel impulse response tracking break- employed to generate branch metncs and a derision (box 

down. By finding the most probable instant at which the 65). This is accomplished usinp the e -gnaTinnQ niit1inp4 

problem might occur, equalizer operation approaches that 65 below. Then it is determ ined if a decision has tvi'-n maHft 

instant from both forward and a time-reversed directions, (box 66). If a decision has been made, then an estimate of the 

both of which begin with known information sequences that channel impulse response is generated (box 68). If th e 
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decision is not made, or once the channe] impulse response 
estimate has been generated, then the symbol number is 
compared to the location of the power fade less a predeter- 
mined number of additional symbols (box 69). Processing is 
then repeated by decrementing the symbol number (box 63) ^ 
and repeating steps (boxes 64-69) until the fade location less 
a predetermined number of additional symbols has been 
reached. 

More particularly, and in operation, samples entering the 
equalizer 21 may be identified as z(n), and the output 
decisions may be identified as a(n). The probability of 
correctness of a(n) depends on location within the bursts. 
When a(n) is known with certainty the values of a(n), 
denoted a,(n), are used by the channel impulse response 15 
estimator 32 for training. At other times, the best estimate of 
a(n) is the output of the traceback decision process , of the 
equalization trellis 31, denoted a^n). 

The equalization trellis 31 operates as follows. Equaliza- 
tion proceeds in the forward direction from the beginning of ^ 
the preamble up until M symbols after the minimum power 
symbol. In the reverse direction, the same occurs with 
processing continuing M symbols beyond the minimum 
power point. This overlap ensures that trace-back through 
the trellis in all likelihood converges to a single path by the 
minimum power point. 

Traceback for actual decisions does not occur until the 
completion of the equalization process. In addition to final 
traceback, however, there is a need for tentative decisions 30 
during equalization, to provide data estimates for the chan- 
nel impulse response estimation to remain current. A trade- 
off in determining these tentative decisions arises (a) 
because the more up-to-date the information is, the more 
up-to-date the chaimel impulse response estimate can be 35 
(remembering that the channel is far from stationary at high 
speeds), and (b) the higher the number of symbols that are 
considered before tentative decisions are made, the more 
accurate the decisions will be; and hence, the lower the 
probability that errors are introduced into the channel ^ 
impulse response estimation. In the case of 4-state equal- 
ization there is very littie sensitivity to the number of 
constraint lengths of delay introduced. 

Branch metrics are calculated in the equalizer 21 using the 
following equation: 



fcrinctric[z(jt), app^tate(J)] = k(k) 



- 2 C,(fc)xaft(/.«)P 
«=0 



Cyjtt^') = Cj,W + Ki [ zik) - ^ Crm X a(* - n) ] a*(Jt - I) 



45 



where app__stateG) represents a hypothetical state in com- 
bination with potential input data; a;,G*ii) is a corresponding 
transmitted signal (constellation point), C represents the 
current estimate of the channel* s impulse response, and z is 
the measured output of the matched filter 27. 

The chaimel estimator 32 utilizes a second order least 
mean square algorithm to determine the coefficients of the 
transversal filter 27 that is an estimate of die channel 
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8 

-continued 



Ci(*+») = Ci W + Cffi(^i> + ^ z{k) - ^ C„(A) X fl(t - n) j - I), 



where Co(k) and Ci(k) arc complex values of estimated 
channel impulse response taps, CsJ^) and C„(k) arc com- 
plex intermediate values related to the estimated channel 
impulse response taps, permitting second order operation, 
Kj and K2 arc the real gain values controlling the tracking 
rate of the channel impulse response estimation process, z(k) 
arc complex symbol spaced sampled outputs of the receiver 
matched filter, and zQs) arc complex estimated or known . 
values of transmitted symbols. 

The values Kj and within these equations control the 
rate of adaptation, and (conversely) the sensitivity to noise 
and decision errors. Consequentiy, to minimize the error 
rate, a trade-off between ability to track changes in the 
chaimel and degradation in performance due to imperfect 
input information is needed to optimize the values of Kj and 
Kj. The optimal values of Kj and Kj vary as a function of 
instantaneous signal to noise ratios, and thus as a function of 
depth of fade. Therefore, algorithms for modifying the 
values during each burst have been evaluated, with consid- 
erable improvement in performance relative to that achiev- 
able with constant settings. 

One approach for modifying Ki and has provided 
good performance and is as follows: 

1. Set the values of Ki and that will apply at die 
symbol determined to correspond to the deepest fade; 
Ki_fade. 

2. Adjust each value linearly (with preset slope — ^Kl_ 
slope and K2_slope) to reach the selected values at Uie fade 
location, using: 

before forward processing — ^initialize 
Ki=Ki_fade-— Ki_slope.fade_location 
^2=^2— fade — Kz^slope.fade — ^location 
before reverse processing — initialize 
Ki=Ki_fade — Ki_slope.( 1 77— fade_location) 
Kj^Kj^fade— K;^slope.(177— fade„location) 
during processing — as each symbol is processed 
K,=Ki+Ki_slope 
K2=K2+K2_slope, 
where Ki_fade is the real value of Ki at the symbol with the 
maximum estimated fade depth, K2_fade is the real value of 
K2 at the symbol with the maximum estimated fade depth, 
Ki_slope is the real increment in K, applied during pro- 
cessing of each symbol, Ka^slope is the real increment in 
K2 apphed during processing of each symbol, and fade_ 
location is the symbol number at the maximum estimated 
fade depth, and last__location is the symbol number of the 
final symbol. 

Estimation of the location of the power fade entails use of 
the received symbols &om the matched filter 27, and the 
settings on the AGC circuit 35 that were active during 
reception of those symbols. As the response of the amplifier 
22 to the AGC circuit settings is effectively instantaneous, 
the primary delays in utilizing this information arise in die 
matched filter 27. This filter 27 is a linear phase filter 
(constant delay), so tiie available input information can be 
easily transformed into an accurate estimate of the envelope 
power This envelope is averaged by a rectangular FIR filter 
over about ten (10) symbol times, with very good perfor- 
mance. 

After completion of acquisition, the carder frequency 
offset should be less than 200 Hz. To operate without 
impairment, diis offset should be on the order of 20 Hz or 
less. Thus, estimation of and correction for carrier offset 
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must continue after acquisition. The method employed uti- 
lizes the fact that when frequency offset occurs, the taps of 
the channel impulse response will rotate consistently at a 
rate proportional to the offset. Changes in tap phases over 
fixed periods, therefore, provide an observable characteristic 
to apply to frequency control. Note that random phase 
changes occur in addition to these consistent rates of change, 
so filtering is used to extract the frequency offset. In practice, 
offsets of around 1000 Hz can be resolved although the 
maximum expected offset after acquisition is 200 Hz. The 
approach used is as follows: 

1 . During the reception of each burst, the haft of that burst 
that does not include the deepest fade is selected for track- 
ing. This scheme is aimed at avoidance of the very high rates 
of change in phase that typically accompany transitions 
through low signal amplitudes. 

2. Two samples of each of the two estimated channel 
impulse response taps are recorded: just after the preamble 
(or leading into the postamble if the fade occurred during the 
first half of the slot), and 20 symbols later (or 20 symbols 
cariier). At a symbol rate of 24300 symbols per second, a 
100 Hz offset would result in an average rotation of 29.6 
degrees during the 20 synibol period. For any rotation in 
excess of 180 degrees, the observed rotation would be less 
than 180 degrees but in the opposite direction. This aliasing 
could impact performance for frequency offsets above about 
300 Hz. In typical operation, however, the detriment to 
performance resulting from such aliasing has proved mini- 
mal, due to the anti-aliasing filtering inherent in the tracking. 
The selection of a sampling window of 20 symbols was 
based on concern about this aliasing. Otherwise, a longer 
window would improve noise immunity. 

3. From information determined during the bit timing fine 
tuning, the dominant tap is selected. Using the recorded 
settings for this tap, a phase change is calculated, yielding an 
estimate of the frequency offset. 

4. These estimates are then filtered over many bursts to 
reduce the "noise" that arises primarily due to the random 
(zero mean) presence of Doppler offsets and Gaussian noise. 
The filter output provides an estimate of the carrier offset 40 
and can be used to directly update the frequency control 
hardware. The offset is given by: 

f_offsel_estiinatejtti=(l-K^^)f_offset_estimatejt+K^^ 
freq_observed, 

where freq_observed is derived from the observed phase 45 
change, the constant K^^ controls the convergence rate of the 
estimation process, f_offset_estimate;t is the estimated fre- 
quency offset at frame "k**, and K^^ is a constant controlling 
the convergence rate of the frequency tracking. If f_offset_ 
estimate reaches half the resolution of the frequency source, 
then a step in frequency is applied, e.g., if the resolution is 
20 Hz and f_offset_estimate exceeds 10 Hz, then a 20 Hz 
change in reference is applied. At the same time f_.offsei_ 
estimate is reinitialize! 

Referring to FIG. 5, it illustrates a flow diagram showing 
the processing performed by the equalizer 20 to iniplement 
carrier frequency offset compensation. Utilizing an already 
located fade, a decision (box 100) is made as to whether to 
use the first or second half of the received slot for firequency 
offset estimation. Based on this decision, samples are taken 
twenty symbols apart in the appropriate half of the slot 
(boxes 101,102). For the seleaed case, individual taps are 
compared and the larger is chosen (decisions 103, 104). The 
phases of the chosen lap at the selected two times are then 
subtracted (boxes 105-108) to produce "fireq_observed", a 
noisy estimate of the offset. Tbiz is filtered (box 109) to 
generate an accurate estimate of the offset If an adjustment 
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in setting of the frequency control would reduce this offset, 
then a decision is made to do so (decision 110); and the 
decision is then implemented (box 111), 

The equalizer is reasonably insensitive to errors in bit 
timing. However, for the following reasons, symbol timing 
adjustments continue during equalizer operation. The initial 
estimate produced by acquisition may differ sufficiently 
from optimal timing so that performance would benefit from 
adjustment. The transmit and receive symbol timing clocks 
may differ by about 5 ppm, resulting in drift of about 0. 1 nS 
per frame (or a symbol every 8 seconds). This drift must be 
compensated for. In practice, individual independentiy-de- 
layed signal paths, will randomly rise and diminish in 
average strength, resulting in situations that would be best 
catered for by different symbol timing. Optimal symbol 
timing depends on an ability to track these changing situa- 
tions. 

The operation of the symbol timing control is as follows! 
The approach has similarities to the early-late gating 
schemes frequentiy employed in direct-sequence spread 
spectrum receivers. As each burst is received, a measure of 
the enor between the expected preamble and the actual 
received preamble is generated. In addition, in alternating 
frames, similar measures are made on time advanced and 
retarded venions of the same input samples. If no timing 
adjustment is necessary, the error generated with the existing 
timing should be less (on average) than eitiier of the others. 
Adjustments are made when this is not the case or there is 
a consistent disparity between the advanced and retarded 
error estimates. This process is simply a search for bit timing 
thai minimizes the error statistic, as iUustrated in FIG. 4. The 
control loop used includes an estimator of any consistent 
change in timing, corresponding to drift with respect to the 
transmitter. Drift in the order of 10 ppm can be compensated 
for by this loop. 

This search for a minimum may be hampered by the 
possible presence of a local (non-global) minimum. In fact, 
for this statistic the presence of two minima is common 
(corresponding to the two taps implicit in the equalizer 
stmcture— see FIG. 1). The approach taken to resolve tiiis 
conflict is as follows. The more advanced minimum is 
presumed to be the preferred sampling time. Multiple 
minima typically arise when there is a small level of delay 
spread, i.e., less than about 10 jiS. Under such conditions the 
ratio of magnitudes of the estimated paths in the (symbol- 
spaced) channel impulse response differs significanUy in the 
region of the more advanced minimum from that in the more 
retarded case. Thus, the ratio of tap magnitudes provides a 
statistic from which to conclude the appropriateness of a 
selected minimum. 

Wth reference to FIGS. 6A and 6B they show flow 
diagrams illustrating the processing performed by the equal- 
izer 20 to implement bit timing control Inputs (box 80) 
include the on-time and time-offset samples (z (n) and z 
offset(n)), and a flag to indicate the direction of the time 
offset. The on-time samples are fed into the equalizer 20 just 
as they are during normal training 83. Similarly, the time 
offset samples are fed to the equalizer 20 (box 84). In both 
cases, the branch metrics (on Uie known correct paths) are 
accumulated over the latter symbols to provide measures 
(ERROR^^ and ERROR OFFSET^ J of the degree to 
which the samples match expectations. 

In a separate process the magnitudes of each of two taps 
estimated as the channel impulse response at the end of the 
traiiung process are calculated (box 85). Averaging the ratio 
of these taps over a immber of frames (boxes 86-89) permits 
a judgement to be made as to whether the bit timing has 



06/08/2004, EAST Version: 1.4.1 



5.537,419 



11 



12 



selected an inappropriate local minimum. If a threshold (box 
90) is reached, then bit liming will be advanced by a ftill 
symbol time (box 91). Taking account of the relative time at 
which samples were taken (box 92), the ERRORc„,„ and 
ERROR OFFSETc„;„ measures are combined to generate a 5 
noisy estimate of an appropriate timing adjustment (boxes 
93, 94). This estimate is then filtered (box 95) to generate an 
actual timing offset adjustment. To compensate for consis- 
tent drift, an additional term "drift_est" monitors and com- 
pensates for this effect. 

Thus there has been described a maximum likelihood 
sequence estimation based equalization method for use in 
mobile digital cellular receivers. It is to be understood that 
the above-described embodiments are merely illustrative of 
some of the many specific embodiments which represent 
applications of the principles of the present inventioa 
Clearly, numerous and other arrangements can be readily 
devised by those skilled in the art without departing from the 
scope of the invention. 

What is claimed is: 

1. In a digital communications receiver in which data is ^ 
received sequentially in lime slots on a transmission chan- 
nel, a method for adjusting the timing of the sampling of the 
data on the transmission channel comprising the steps of; 

sampling the data received on the transmission channel ^ 
sequentially with a determined adjustable timing; 

storing the samples taken of the transmission channel 
data; 

generating from the stored samples, a first estimate of a 
transmission channel impulse response and a first error 30 
measurement comprising a measure of the degree to 
which the first transmission channel impulse response 
estimate combined with known data varies from the 
stored samples; 

generating, from the stored samples, a second estimate of 
a transmission channel impulse response and a second 
error measurement comprising a measure of the degree 
to which the second transmission channel impulse 
response estimate combined with the known data varies 
from the stored samples, the second error measurement 40 
using timing that is offset relative to the timing for the 
estimate for the first transmission channel impulse 
response estimate; and 

determining, by combining the first, and second error 
measurements, a sample timing that minimizes the 
combined error measurement and adjusting the sample 
timing in accordance with the determination. 

2. The method of claim 1 wherein the received data 
corresponds to the known data transmitted into the trans- 
mission channel to the receiver and wherein the channel 
impulse response estimates are generated using a maximum 
likelihood sequence estimation process, the maximum like- 
lihood sequence estimation process using variable taps that 
are determined by: 

generating a set of synthesized samples by applying a 
portion of the known data corresponding to the stored 
samples of the received data to an estimated channel 
impulse response; 

determining the square of the difference between the 50 
stored samples and the set of synthesized samples; 

modifying the taps of the estimated channel impulse 
response and generating a second set of synthesized 
samples by applying the portion of the known data 
conesponding to the stored samples of the received 65 
data to the modified estimated channel impulse 
response; 
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repeating the steps of modifying the taps and generating 
a set of synthesized samples and determining the square 
of the difference in an iterative manner in order to 
minimize the square of the difference between the 
stored samples and the synthesized samples. 

3. The method of claim 2 fimher comprising the steps of: 
determining niagnitudes of each of the taps; 
averaging a ratio of the magnitudes of the taps over a 

prescribed number of frames; 

determining whether the ratio of the magnimdes of the 
taps having been averaged over the prescribed number 
of frames has exceeded a prescribed threshold; and 

advancing, in the event the ratio of the magnimdes of the 
taps having been averaged over the prescribed number 
of frames has exceeded the prescribed threshold, the 
sample timing by a prescribed amount of time. 

4. In a digital conMnunications receiver in which data is 
received sequentially in time slots on a transmission chan- 
nel, a sample timing adjuster for the timing of the sampling 
of the data on the transmission channel comprising; 

an analog to digital converter for sampling the data 
received on the transmission channel sequentially with 
a determined adjustable timing; 

a memory for storing the samples taken of the transmis- 
sion channel data; 

means coupled to the memory for generating, from the 
stored samples, a first estimate of a transmission chan- 
nel impulse response and a first error measurement 
comprising a measure of the degree to which the first 
transmission channel impulse response estimate com- 
bined with known data varies from the stored samples; 

means coupled to the memory for generating, from the 
stored samples, a second estimate of a transmission 
channel impulse response and a second error measure- 
ment comprising a measure of the degree to which the 
second transmission channel impulse response estimate 
combined with the known data varies from the stored 
samples, the second error measurement using timing 
that is offset relative to the timing for the estimate for 
the first transmission channel impulse response esti- 
mate; and 

means coupled to the first, and second estimate generating 
means and to the analog to digital converter for deter- 
niining, by combining the first, and second error mea- 
surements, a sample timing that minimizes the com- , 
bined error measurement, and for adjusting the sample 
timing in accordance with the determination. 

5. The adjuster of claim 4 wherein the received data 
corresponds to the known data transmitted into the trans- 
mission channel to the receiver and wherein each of the 
means for generating charmel impulse response estimates 
comprise means for using variable taps as part of a maxi- 
mum likelihood sequence estimation process, said means for 
using variable taps comprising means for performing the 
following functions: 

generating a set of synthesized samples by applying a 
portion of the known data corresponding to the stored 
samples of the received data to an estimated channel 
impulse response; 

determining the square of the difference between the 
stored samples and the set of synthesized samples; 

modifying the taps of the estimated channel impulse 
response and generating a second set of synthesized 
samples by applying the portion of the known data 
corresponding to the stored samples of the received 
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data to the modified estimated channel impulse 
response; 

repeatedly modifying the taps and generating a set of 
synthesized samples and determining the square of the 
difference in an iterative manner in order to minimize 5 
the square of the difference between the stored samples 
and the synthesized samples. 
6. The adjuster of claim 5 wherein the received data 
coiresponds to the known data transmitted into the trans- 
mission channel to the receiver and wherein each of the 10 
means for generating channel impulse response estimates 
comprise means for using variable taps as part of a maxi- 
mum likelihood sequence estimation process, said means for 
using variable taps further comprising means for performing 
the following functions: 
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determining magnitudes of each of the taps; 

averaging a ratio of the magnitudes of the taps over a 
prescribed number of frames; 

determining whether the ratio of the magnitudes of the 
taps having been averaged over the prescribed number 
of frames has exceeded a prescribed threshold; and 

advancing, in the event the ratio of the magnitudes of the 
taps having been averaged over the prescribed number 
of frames has exceeded the prescribed threshold, the 
sample timing by a prescribed amount of time. 

♦ ♦ * ♦ * 
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